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Abstract: The mass of the strange quark is determined from SU(3)-breaking effects 
in the r hadronic width. Compared to previous analyses, the contributions from scalar 
'-^ and pseudoscalar spectral functions, which suffer from large perturbative corrections, are 

replaced by phenomenological parametrisations. This leads to a sizeable reduction of the 
uncertainties in the strange mass from r decays. Nevertheless, the resulting value is 
5^ , still rather sensitive to the moment of the invariant mass distribution which is used for 
the determination, as well as the size of the quark- mixing matrix element \Vus\- Imposing 
the unitarity fit for the CKM matrix, we obtain ms{2 GeV) = 117 ± 17MeV, whereas for 
the present Particle Data Group average for \Vus\, we find ms(2GeV) = 103 ± ITMeV. 
On the other hand, using an average of from other sources as an input, we are able to 
calculate the quark-mixing matrix element \Vus\, and we demonstrate that if the present 
measurement of the hadronic decay of the r into strange particles is improved by a factor 
of two, the determination of \Vus\ is more precise than the current world average. 
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1 Introduction 



About a decade after the development of QCD sum rules by Shifman, Vainshtein and 
Zakharov it was realised that the hadronic decay of the r lepton could serve as an 
ideal system to study low-energy QCD under rather clean conditions ^. Shortly 

afterwards, the experimental precision of r decay data was sufficiently improved in order 
to explore this possibility in practice. Since then, detailed investigations of the r hadronic 
width as well as invariant mass distributions have served to determine the QCD coupling 
as to a precision competitive with the current world average More recently, the 

experimental separation of the Cabibbo-allowed decays and Cabibbo-suppressed modes 
into strange particles opened a means to also determine the mass of the strange quark 
l^,|TI],[l^,|l3|,[l^,|l5|,|l^,|l3, one of the fundamental QCD parameters within the Standard 



Model. 

Until today, these strange mass determinations suffer from sizeable uncertainties due to 
higher order perturbative corrections. In the sum rule under investigation, scalar and pseu- 
doscalar correlation functions contribute, which are known to be inflicted with large higher 
order QCD corrections 0, [T8|, |19|, [1^] , and these corrections are additionally amplifled by 
the particular weight functions which appear in the r sum rule. As a natural continuation, 
it was realised that one remedy of the problem would be to replace the QCD expres- 
sions of scalar and pseudoscalar correlators by corresponding phenomenological hadronic 
parametrisations [l^, ^ , which are expected to be more precise than their QCD 

counterparts. In this work, we shall present a complete analysis of this approach, and it 
will be shown that the determination of the strange quark mass can indeed be signiflcantly 
improved. 

By far the dominant contributions to the pseudoscalar correlators come from the kaon 
and the pion, which are very well known. The corresponding parameters for the next 
two higher excited states have been recently estimated |^. Though much less precise, 
the corresponding contributions to the r sum rule are suppressed, and thus under good 
theoretical control. The remaining strangeness-changing scalar spectral function has been 
extracted very recently from a study of S-wave Kit scattering [^, ^] in the framework 
of chiral perturbation theory (xPT) [^,|2^ with explicit inclusion of resonances [^6|,^ . 



The resulting scalar spectral function was then employed to directly determine the strange 
quark mass from a purely scalar QCD sum rule ||28|. On the other hand, now we are also 
in a position to incorporate this contribution into the r sum rule. The scalar ud spectral 
function is still only very poorly determined phenomenologically, but it is well suppressed 
by the small factor (m„ — m^)^ and can be safely neglected. 
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In the next section, for completeness, we begin by recalling the theoretical expressions 
which are required in the sum rules under investigation. Then, in section 3, the scalar and 
pseudoscalar contributions to the r sum rules are investigated separately, and it is demon- 
strated that, although the uncertainties on the QCD side are large, they are well satisfied. 
In section 4, an improved determination of the mass of the strange quark is performed by 
replacing the QCD expressions for the scalar and pseudoscalar spectral functions with the 
corresponding hadronic parametrisations. Since the resulting strange quark mass depends 
quite sensitively on the value of the quark- mixing matrix element \Vus\, in section 5 we 
take a different approach to the considered sum rule. 

Assuming a strange mass as extracted from other recent analyses, we use the r sum 
rule to determine \Vus\- The resulting uncertainty on \Vus\ is still sizeable, but the error 
is dominated by the partial hadronic width of the r decaying into strange particles. It 
is then demonstrated that an improvement of this measurement by a factor of two would 
result in a determination of \ Vus\ better than the present world average, thus shedding some 
light on the question of the unitarity violation in the first row of the Cabibbo-Kobayashi- 
Maskawa (CKM) matrix. Finally, in section 6 we compare our results with other recent 
determinations of the strange quark mass and in section 7 we end with a discussion and 
conclusions. 



2 Theoretical framework 

Below, we review the main theoretical expressions required in the theoretical analysis of 
the inclusive hadronic r decay width. Further details and complete expressions can be 



found in the original works p,[Ty,[T^. The central quantities in such an analysis are the 
two-point correlation functions 

KiUp) = ' I dxe'^^ {n\ T{J^ifix) (0)t}|fi) , (2.1) 

where fl denotes the physical vacuum and the hadronic vector /axialvector currents are 
given by J^^{f{x) = {(ijlfi{l5)qi){x). The indices i,j denote the light quark fiavours up, 
down and strange. The correlators Il^l/^j{p) have the Lorentz decomposition 

KIUp) = iP^P^ - 9,up') Tlll^^^p') + P,Pu n.y^'^p^) , (2.2) 



where the superscripts in the transversal and longitudinal components denote the corre- 
sponding angular momentum J = 1 (T) and J = (L) in the hadronic rest frame. 
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The hadronic decay rate of the r lepton, 

_ r(r~ ^ hadrons i^r(7)) 

Rr = 



Rt,V + Rt,A + Rt,S 



(2.3) 



r(r ^ e-iyeJ^rh)) 

can be expressed as an integral of the spectral functions Imn-^(s) and Imn^(s) over the 



invariant mass s = of the final state hadrons ||29|| : 



127V 



' as / s Y 



M2 



1 + 2— )lmn^(s) + Imn^ 

Mi 



The appropriate combinations of the two-point correlation functions are given by 



IK 



ud 



n 



V,J 
ud ' 



s) + n; 



A,J, 
ud \ 



AJ, 



(2.4) 



(2.5) 



with Vij being the corresponding matrix elements of the Cabibbo-Kobayashi-Maskawa 
quark- mixing matrix. As has been indicated in eq. ( ^.31 ), experimentally, one can disen- 
tangle vector from axialvector contributions in the Cabibbo-allowed {ud) sector, whereas 
such a separation is problematic in the Cabibbo-suppressed {us) sector, since G-parity is 
not a good quantum number in modes with strange particles. 

Additional information can be inferred from the measured invariant mass distribution 
of the final state hadrons. The corresponding moments R^\ defined by 

r T 

can be calculated theoretically in analogy to Rr = R^ . Exploiting the analytic properties 
of n'^(s), the moments ( p. 61) can be expressed as contour integrals in the complex s-plane 

2. 



running counter clockwise around the circle |s 

dx r 

X 

\x\=l 



R. 



kl 



ITT 



^J^L+A^) V'^^' {M'^x) +AJ^1\x) P^(M^a; 



L+T 



kl/ 



(2.7) 



Here, we have defined the new integration variable x = s/M^, and general expressions 
for the kinematical kernels ^^^^(x) and J^l{x) have been presented in ref. [0. For the 
convenience of the reader, we have compiled explicit expressions for the moments (0, 0) to 
(4, 0), which will be utilised in our phenomenological analysis in table 

In the derivation of the above expression (|2.7| ), integration by parts has been employed 
to rewrite R^'' in terms of the physical correlators V^^^{s) and T>^{s), 



V 



L+T, 



S) = —s 



d_ 

ds 



n 



L+T / 



/ N s d 



(2.8) 



3 



{k,l) 




^f(x) 


(0,0) 
(1,0) 
(2,0) 
(3,0) 
(4,0) 


{l-xf{l + x) 
^(1 -x)4(7 + 8x) 
^ (1 - x f (4 + 5x) 
i (1 - x)6 (3 + 4x) 

n(i-a;)M5 + 7x) 


(1 -a;)3 
f(l-x)^ 
f (1-^f 

1(1 -^r 

f (l-x)7 



Table 1: Explicit expressions for the kinematical kernels J^l^tI^) J-'l{x) which will 
be used in our phenomenological analysis. 



which both satisfy homogeneous renormalisation group equations, and thus eliminate the 
dependence on unphysical (renormalisation scale and scheme dependent) subtraction con- 
stants. For large enough negative s, the contributions to V-^{s) can be organised in the 
framework of the operator product expansion (OPE) in a series of local gauge-invariant 
operators of increasing dimension D = 2n times appropriate inverse powers of s. This 
expansion is expected to be well behaved along the complex contour \s\ = M^, except 
at the crossing point with the positive real axis As can be seen from eq. (|2.7|) and 
table |l|, however, the contribution near the physical cut at s = is strongly suppressed 
by a zero of order three or larger. Therefore, uncertainties associated with the use of the 
OPE near the time-like axis are expected to be very small. 

Inserting the OPE series for {s) into ( p.7|) and performing the contour integration. 



the resulting expression for Elf can be written as 
El 



= 3 (|K.r + \Vus\')SeJ 1 + S'^^'^ + E (cos^^c^^r + sm'ec6T^)] (2.9) 

D>2 ^ 

where sin^^c = |KjsP/(|V^dP + iKisP), (^c being the Cabibbo angle, and the electroweak 
radiative correction Sew = 1.0201 ±0.0003 [p2| , |3^ , |3^ has been pulled out explicitly. Owing 
to chiral symmetry, the purely perturbative dimension-zero contribution S''''^^^ is identical 
for vector and axialvector parts. The symbols 6^^^^^ = i^udy + ^uIIa)/'^ stand for the 
average of vector and axialvector contributions from dimension D > 2 operators which 
contain imphcit suppression factors 1/M^. 

The separate measurement of Cabibbo-allowed and Cabibbo-suppressed decay widths 



of the r lepton |12] allows one to pin down the flavour SU(3)-breaking effects, dominantly 
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induced by the strange quark mass. Defining the difference 

5R- ^ %^ - ^ = 3 ^Ew E (^T^ - ST') , (2.10) 

many theoretical uncertainties drop out since these observables vanish in the SU(3) hmit. 
In particular, they are free of possible flavour-independent instanton as well as renormalon 
contributions which could mimic dimension-two corrections. 

As was already stated in the introduction, the longitudinal contributions S'^j^l' are 
plagued with huge perturbative higher order corrections, which in previous analyses re- 
sulted in large corresponding uncertainties for the strange quark mass. Thus, in the follow- 
ing, we shall replace the theoretical expressions for the longitudinal part by corresponding 
phenomenological contributions, thereby strongly reducing those uncertainties in the 
determination. Before we embark on the full investigation of the moment differences 6R'^'', 
to be able to judge the gain of the intended replacements in the longitudinal sector, in the 
next section we shall flrst investigate r-like sum rules for the longitudinal part alone. 



3 The longitudinal sector 

From eq. (|2.7|) , again making use of the analytic structure of the correlation functions, one 
can deduce the following r-like sum rules for the longitudinal contributions: 

^S#M = -24vr^/rfx(l-x)^+V+Vr''(^'^) = -4vr^/ ^ (x) D^^'^M^) , 

|x|=l 

(3.1) 

with Pij^'^{s) = Imn|^^'^'^(s)/7r being the longitudinal spectral functions. Note that for 
convenience, we have defined Rijy/j^ without the appropriate CKM and electroweak correc- 
tion factors. The longitudinal correlators are directly related to corresponding correlators 
^J^^'^(s) for the divergences of vector and axialvector currents: 



nr''(^) = ^[^r(^)-^r(o)j- (3-2) 
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Furthermore, via the equations of motion for the quark fields, ^|^^^(s) is related to the 
two-point correlators 11^^ ^(s) for the scalar /pseudoscalar currents jf-^{x) = {qj{'^^)qi){x) , 
and the constants \l/]^'''^(0) to the quark condensate from a Ward-identity |3 

vl/J/^(s) = {m,Tm^fli^l\s), vI;J/^(0) = - (m, ^ m,)[ ^ ] , (3.3) 



5 



where and rrij are the masses of the quark flavours qi and qj respectively, and {(jiqi) as 
well as (qjqj) are the corresponding quark condensates. The upper sign corresponds to the 
vector/scalar and the lower sign to the axialvector/pseudoscalar case. For simplicity, the 
vacuum state has been omitted in the condensates, and we would like to point out that the 
bi-quark operators in eq. ( |3.3|) are non-normal-ordered |^,^. The theoretical expressions 



for the longitudinal correlators '^'^^•(s) can either be obtained from refs. P,|T8|,|T0|,|T3|, or can 
be calculated straightforwardly from the second of eqs. (|2.^) together with eq. (|3.2|) and the 



expressions for \E']^''^(s) given in ref. |3^. For further details and the explicit expressions 



the reader is referred to these references. 



Theory: 
Phenom: 



00,L 

us, A 



O0,L 
us,V 



pOO,L 
''^ud,A 



-0.144 ±0.024 
-0.135 ±0.003 



-0.028 ±0.021 
-0.028 ±0.004 



-(7.79 ±0.14) ■ 10-3 
-(7.77 ±0.08) • 10-3 



Table 2: Comparison of theoretical and phenomenological longitudinal contributions to the 
(0,0) moment of the r sum rule. 



Let us now proceed with a numerical investigation of the sum rules of eq. ( |3.1| ). In 
table 0, we display a comparison of theoretical and phenomenological values of the longi- 
tudinal contributions -R°°'^, R^g'v -^ndA rule. Due to the global mass 
squared factor, these results depend sensitively on the quark masses. For definiteness, 
we have set the strange mass to ms(2GeV) = 105 MeV, a value compatible with most 

^,|28| and lat- 



recent determinations of both from QCD sum rules p9|, |40|, p!5|, p!4|, iTi 
tice QCD |^l|,^,^^. A detailed comparison to other rris determinations will be pre- 
sented in section 6. The light masses and have been fixed from the xPT ratios 
^u/^d = 0.553 ± 0.043 and rris/md = 18.9 ± 0.8 |^5[. On the theoretical side, we have 
also included instanton contributions according to the model used in , because they are 
known to play some role in scalar /pseudoscalar finite energy sum rules. 

The theoretical uncertainties in table ^ are largely dominated by higher order pertur- 
bative corrections. Already for the (0, 0) moment, the first and second order corrections 
are of similar size, while the third order correction is larger than these two. Thus for 
this moment, in the spirit of asymptotic series, we have cut the perturbative series at the 
second order. For higher moments, the behaviour is still worse and there the second order 
turns out to be larger than the first one, which then entails even larger perturbative un- 
certainties. For this reason, we have not presented detailed results for the higher moments 
although also in these greement between theoretical and phenomenological results 
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is found within the errors. A detailed discussion of how we estimated the theoretical errors 
will be given in the next section. 

The uncertainties in the phenomenological numbers, on the other hand, are much 
smaller. For the pseudoscalar spectral functions we have employed the parametrisation 
by Maltman and Kambor ||2l| . In the case of the pseudoscalar us channel, the correspond- 
ing expression reads 



2/|-M^5(s - Ml) + 2flMtB,{s) + 2flM^B^{s) , 



(3.4) 



with fi and Mj being decay constants and masses of the higher exited resonances. Bi[s) is 
the standard Breit-Wigner resonance shape function 



BAs] 



(3.5) 



For the pseudoscalar ud channel an equation analogous to ( p.4|) holds. The masses and 
widths of the resonances have been taken from the Review of Particle Physics and the 



decay constants from ref . . For convenience, we have compiled these input parameters 
in table ^. The dominant uncertainties in the phenomenological results of table ^ are 
due to the errors in the decay constants fx and /^r, as well as the variation of the higher 
resonance decay constants as given in table |^. Changing resonance masses and widths 
within reasonable ranges only has a minor effect. Therefore, for these quantities, we have 
not given explicit uncertainties. 





7r(1300) 


7r(1800) 


K{lAm) 




Mi [MeV] 


1300 


1800 


1460 


1830 


r. [MeV] 


400 


210 


260 


250 


/. [MeV] 


2.20 ±0.46 


0.19 ±0.19 


21.4 ±2.8 


4.5 ±4.5 



Table 3: Input parameters for the pseudoscalar resonance parametrisation. The masses 
and widths are taken from [^ , whereas the decay constants have been determined in [pi 



The scalar us spectral function has been calculated in ref. [|2^ . Below 2 GeV the 
dominant hadronic systems which contribute in this channel are the Ktt and Krj' states 
with ^^0(1430) being the lowest lying scalar resonance. The scalar us spectral function can 
then be parametrised in terms of the scalar Fk-k{s) and FKn'{s) form factors. These form 
factors were obtained in from a coupled-channel dispersion-relation analysis. The S- 



wave Kit scattering amplitudes which are required in the dispersion relations were available 
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from a description of S-wave Kii scattering data in the framework of unitarised xPT with 
2^ . The dominant uncertainty in the scalar us spectral function is due to an 



resonances 



integration constant which emerges while solving the coupled channel dispersion relations. 
In p3[ it was decided to fix this constant by demanding FKniM]^ — M^) = 1.22 ± 0.01, a 



result which from xPT is known to be very robust, since higher-order chiral corrections are 
well suppressed . The uncertainty in the phenomenological value for R^g'y in table ^ 
thus corresponds to a variation of this parameter within the given range. Further details 



on the scalar us spectral function can be found in ref. ||28|| . 

Generally, from table |^ one observes that the longitudinal sum rules are very well 
satisfied. This is partly due to the fact that the subtraction constant "^udi^) dominates the 
ud sum rule and also in the pseudoscalar us sum rule "^^si^) gives a large contribution. In 
the leading order in xPT, these subtraction constants are fixed by the pion and kaon pole 
contributions which led to the famous Gell-Mann-Oakes-Renner (GMOR) relations p7| . 
Therefore, we have a well known piece which appears on both sides of the sum rule. 
Nevertheless, the agreement within errors is non-trivial and also points to the fact that 
the used light quark mass values are reasonable. However, the complete agreement in the 
scalar us case should be considered as accidental. From table || it becomes obvious that 
the theoretical uncertainties in the flavour-breaking r sum rule can be greatly improved by 
replacing the theoretical longitudinal contributions by the corresponding phenomenological 
ones. 

The observed sensitivity of the longitudinal sum rules on the subtraction constants 
allows to obtain further information on these quantities. For this purpose it is convenient 
to rewrite the subtraction constants as 

^tM = 2flMl ( 1 - 5J , ^t{^) = 2/iMi ( 1 - 5k ) , (3.6) 

where the leading term corresponds to the GMOR relations, and to actually determine the 



higher order chiral corrections and 6k l^liEU- addition, from the scalar subtraction 
constant \E'^g(0), the flavour breaking ratio of the quark condensates (ss) / {uu) can be 
estimated. Now we have to also vary the strange quark mass, and we chose ms(2 GeV) = 
105 ± 20 MeV, the central value already advocated above with an uncertainty which is 
slightly more generous than the ranges obtained in the recent scalar and pseudoscalar sum 
rule determinations [^,^. 



Again varying all other input parameters according to the ranges presented here and 
in the next section, we obtain: 

(ss) 

6^ = 0.049 ±0.030, 6k = 0.50 ±0.23, = 0.80 ±0.58. (3.7) 

{uu) 
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These results are in very good agreement to the findings of ref. |^8|, where the values 
of 5^ = 0.047 ± 0.017 and 5k = 0.61 ± 0.22 were estimated in the framework of xPT 
assuming a given input for the ratio {ss)/{uu). However, the latter results of are 
somewhat more precise, since the uncertainties in eq. (|3.7|) are again dominated by the 
large higher order perturbative corrections and now in addition by the error on the strange 
quark mass. Besides from the scalar subtraction constant, the ratio (ss) / (uu) can also be 



directly inferred from Sk, and from 6-,^ by inverting the line of reasoning applied in ref. pH 



The central values found in this way are {ss)/{uu) = 0.63 and 0.75 respectively, however 
with uncertainties which are even larger than the one presented in eq. (|3.7| ). Nevertheless 
it is gratifying to observe that the central results for {ss) / {uu) display good consistency. 



4 Strange quark mass 

Our determination of the strange quark mass from the SU(3)-breaking differences of 
eq. ( |2.10| ) proceeds in complete analogy to the previous analyses |TB|,|T^ , with the exception 
that now we have replaced the longitudinal contributions to the sum rule by phenomeno- 
logical parametrisations as discussed above. Solving for the strange mass, eq. ( p.lOj ) can 
be written as 



''Tjphen 



18(1 - e^)Afc+ (a,) |^ 5ew 



where = rrid/ms. The phenomenological contribution containing the longitudinal spec- 
tral functions is given by 

r Tjkl^L Tjkl,L . -r)kl,L jjkl^L -Qkl,L /a r\\ 

'^^T,p\ica — ^ud.y + ^ud,A ~ ^us,V ~ ^us,A i l^-^J 

with Rfjv/A defined in eq. (|3.1|) and calculated using the scalar and pseudoscalar spectral 
functions discussed in the last section. 

The perturbative QCD correction A^j^^{ar), associated with the dimension-2 contribu- 
tion, depends only on a,- = as{MT-)/n and has been discussed in great detail in refs. |]10|,|13| . 



Since we have subtracted the badly behaved longitudinal contribution, the remaining part 
displays a better convergence at higher orders. In contour-improved perturbation the- 
ory |^9|, and taking as a central value ^^(Mt-) = 0.334 ||^, the expansion for the relevant 



moments takes the form 



A^o+^(a^) = 0.753 + 0.214 + 0.065 - 0.051 + ... = 0.981 
Afo+^(a^) = 0.912 + 0.334 + 0.192 + 0.056 + ... = 1.493 
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A^+^(a^) = 1. 055 + 0.451 + 0.330 + 0. 189 + ... = 2.024, (4.3) 
^fo^^K) = 1.190 + 0.571 + 0.484 + 0.352 + ... = 2.597, 
^40^ (ar) = 1.324 + 0.697 + 0.657 + 0.551 + ... = 3.228. 

The last term in the series, corresponding to the order correction, has not yet been 
calculated. In eq. (^73|) , we have employed an estimate based on assuming a geometrical 
growth of the Adler function V'^^'^{s) |l^. Finally, 5E^'^^'^ are higher-dimensional op- 



erator corrections in the framework of the operator product expansion, the dominant one 
being the D = A correction resulting from the quark condensate. Explicit expressions for 
these contributions can be found in ref. |]13|. In our numerical analysis below, the Z) = 4 



light up and down quark mass corrections as well as the D = Q condensates have also been 
included, although their influence on the value of is insignificant. 



Parameter 


Value 


(0,0) 


(1,0) 


(2,0) 


(3,0) 


(4,0) 


m,(M,) 








197.3 


164.0 


136.6 


115.2 


98.6 




1 


17 




+52.7 
-67.7 


+22.2 
-25.4 


+13.1 
-14.5 


+9.3 
-10.2 


+7.3 
-7.9 


IK.I 


0.2225 + 0.0021 


+28.8 
-32.6 


+14.6 
-16.0 


+9.4 
-10.1 


+6.7 
-7.1 


+5.3 
-5.6 




2xO(af) 
no 0{al) 


+5.9 
-5.4 


-3.1 
+3.3 


-6.0 

+6.9 


-7.1 
+8.6 


-7.5 
+9.6 




1.5 
0.75 




-11.7 
+29.0 


-1.5 
+11.8 


+3.5 
+4.2 


+6.3 
-0.1 


+7.9 
-2.7 


«,(M,) 


0.334 + 0.022 


+14.2 
-9.8 


+5.0 
-3.0 


+1.1 
+0.1 


-0.9 
+1.8 


-2.0 
+2.8 


{ss) / (uu) 


0.8 


+ 0.2 


-0.8 
+0.8 


-3.6 
+3.5 


-5.6 
+5.3 


-7.2 
+6.8 


-8.4 
+7.9 


„scalar 
rus 


see text 


-2.0 
+1.8 


-0.8 
+0.7 


-0.4 
+0.3 


-0.2 
+0.2 


-0.1 
+0.1 


Ik 


113 + lMeV 


-1.1 
+1.1 


-0.8 
+0.8 


-0.6 
+0.6 


-0.5 
+0.5 


-0.5 
+0.4 


fK{iAm) 


21.4 + 2. 


8MeV 


-1.1 
+0.9 


-0.4 
+0.4 


-0.2 
+0.2 


-0.1 
+0.1 


-0.1 
+0.1 


Total 








+68.4 
-76.9 


+29.9 
-30.6 


+18.9 
-19.5 


+17.2 
-16.0 


+17.5 
-15.2 



Table 4: Central results for ms{Mr) corresponding to the unitarity fit for \Vus\ extracted 
from the different moments, as well as ranges for the main input parameters and resulting 
uncertainties for nis- For a detailed description see the discussion in the text. 



Due to cancellations in the difference (|2.10|) , it was already found in the earlier analyses 
l^,|l^|2l| that the strange mass resulting from eq. (|4.1j ) depends sensitively on the value 
for \Vus\- Therefore, in what follows we shall perform two separate extractions of m^. In 

the value for \Vus\ is taken from the Particle Data Group 



the first one, like in ref. 17 



unitarity fit of the CKM matrix [||] yielding = 0.2225 + 0.0021. In table |, we show 
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a detailed account of our results. The first row displays the values of ms{Mr) obtained 
from the different moment sum rules and central values for all input parameters. In the 
following rows, we have listed those input parameters which dominantly contribute to the 
uncertainty on m^, the ranges for these parameters used in our analysis and the resulting 
shift in irts. Only those parameters have been included in the list which at least for one 
moment yield a shift of larger than 1 MeV. Finally, in the last row, we display the total 
error that results from adding the individual uncertainties in quadrature. 

Let us now discuss the different uncertainties in more detail. The experimental results 
for 5B}!f with k = 0...4, as well as the errors corresponding to a variation of \Vus\, 
have been presented in ref. and are taken over from there. They dominate for low 
k and become smaller if k is increased. For the uncertainties related to higher order 
perturbative corrections we have considered three different sources. As already mentioned 
above, the order correction to A^^^{ar) has not yet been calculated and thus we rely 
on estimates of this correction [§,|4^. To be conservative, for our error estimate we have 
chosen two extreme cases, one where this correction is doubled and the other where it 
is removed completely. The corresponding variations of are given in table ^. The 
second source is the renormalisation scale dependence which arises due to the missing 
higher orders which would make physical quantities independent of the scale. Like in 
refs. [^, we have parametrised this dependence by the parameter ^ = /i/a/^, and we 
have varied this parameter in the range 0.75 < ^ < L5. The upper range is slightly 
reduced compared to the previous analyses, but since we have included exphcitly the 
variation of the 0{al) correction, we consider this reduction justified. Finally, our input 
value for ^^(Mt-) = 0.334 ± 0.022 [T^ reflects itself in an additional uncertainty for the 
determination. 

The only parameter which matters for the uncertainty related with the D = 4 contri- 
bution in the OPE is the flavour-breaking ratio of the quark condensates {ss)/{uu). For 
this parameter we have chosen the range (ss)/ (uu) = 0.8 ± 0.2 which includes most values 
found in the literature ^,[5^,^,0. Certainly, this value is also compatible with the 
results obtained in section 3. For increasing k, higher order terms in the OPE become more 
important, and thus for the larger k values this error, together with the perturbative QCD 
uncertainty, dominates. Although we have included the D = 6 contribution according to 
the expressions given in |jl3l in our numerical analysis, for the resulting value of and 
thus also for the uncertainty they can be neglected. 

The final uncertainties included in table ^ concern the variation of the phenomenological 
scalar and pseudoscalar spectral functions. The variations of the corresponding parameters 
have already been discussed in section 3, where we have analysed r-like sum rules for the 
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longitudinal contributions separately. The dominant sources of error are the scalar us 



spectral function as discussed briefly in section 3, and in more detail in ref. as well as 
the decay constants in the pseudoscalar us channel fx and /i^(1460). However, as can be 
seen from table ^, the influence on the error of the strange quark mass is very minor. In 
addition, it again corroborates the reduction of the final uncertainty of achieved through 
the replacement of theoretical by phenomenological scalar and pseudoscalar longitudinal 
contributions. Variations of all remaining input parameters lead to changes in of less 
than 1 MeV, and have thus been neglected. 

To calculate a final result for ms(M^), we have used a weighted average over the different 
moments where as the individual error for each moment, we have taken the larger one. This 
prescription then yields: 

ms{Mr) = 121.7 ± 17.2 MeV => m,(2GeV) = 117 ±17 MeV. (4.4) 

The result on the right-hand side is the value of rris evolved to a renormalisation scale of 
2 GeV. As the uncertainty we have taken the smallest error of one individual moment, 
(3,0) in this case. Taking the correlation of the different moments into account, the final 
uncertainty could still be slightly reduced. However, since the moments are strongly cor- 
related, this would give little improvement, and thus at present we have not incorporated 
this approach. 

Comparing our central strange mass average to the rris values obtained for the individual 
moments, it is evident that there is a strong dependence of rris on the number of the moment 



k. For our approach, this dependence is even stronger than in the analysis of ref. [|T^. The 
reason lies in the fact that the scalar and pseudoscalar sum rules by themselves are satisfied 
for smaller strange masses than the full r sum rule. Thus, when replacing the theoretical 
scalar and pseudoscalar spectral functions by phenomenological ones, the resulting 
moves up. Nevertheless, one should clearly state that the experimental uncertainties, 
especially for the low moments are still large, and the difference between the individual 
moment results and the average strange mass for all values lies in the range of one standard 
deviation. 

As was already mentioned above, the central PDG average for = 0.2196 ± 0.0026 
which is based on the analyses lies more than one sigma below the value 



extracted by requiring unitarity for the first row of the CKM matrix. Due to this large 
difference, we consider it justified to repeat our strange mass determination for this value 
of \Vus\ separately. To properly perform the respective analysis, one would in principle 
need the experimental results for the i?-ratios R^^yj^j^ and R^^g for all moments separately. 
However, as yet, these are only published for the lowest (0, 0) moment, whereas for the 
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higher moments the combined result 5R^ based on the unitarity value for \Vus\ is solely 
available |17|. Since in ref. |17] the uncertainty from varying IV^^I = 0.2225 ± 0.0021 
was listed separately, we can nevertheless work around this problem by rescaling 5R^^ 
appropriately to the above PDG average. 



Parameter 


Value 


(0,0) 


(1,0) 


(2,0) 


(3,0) 


(4,0) 






150.8 


141.5 


122.5 


105.3 


90.8 




experiment 


+61.9 
-103.3 


+25.0 
-30.1 


+14.5 
-16.4 


+10.2 
-11.3 


+7.9 
-8.7 


IK.I 


0.2196 ±0.0026 


+42.1 
-56.6 


+20.3 
-23.5 


+12.7 
-14.2 


+8.9 
-9.8 


+7.0 
-7.6 




2xO(a3) 
no e>(a3) 


+4.3 
-3.9 


-2.6 
+2.8 


-5.3 
+6.1 


-6.5 
+7.9 


-6.9 
+8.9 


e 


1.5 
0.75 


-8.1 
+20.2 


-1.1 
+9.7 


+3.2 
+3.6 


+5.8 
-0.1 


+7.3 
-2.5 




0.334 ±0.022 


+10.1 
-7.0 


+4.2 
-2.5 


+1.0 
+0.1 


-0.8 
+1.6 


-1.9 
+2.5 


{ss) / {mi) 


0.8 ±0.2 


-1.2 
+1.1 


-4.3 
+4.1 


-6.3 
+5.9 


-8.0 
+7.4 


-9.2 
+8.5 


„scalar 
rus 


see text 


-2.5 
+2.2 


-0.9 
+0.8 


-0.4 
+0.4 


-0.2 
+0.2 


-0.2 
+0.1 


h 


113±lMeV 


-1.4 
+1.4 


-0.9 
+0.9 


-0.7 
+0.7 


-0.6 
+0.6 


-0.5 
+0.5 


/ft:{1460) 


21.4 ± 2.8 MeV 


-1.3 
+1.2 


-0.5 
+0.4 


-0.2 
+0.2 


-0.1 
+0.1 


-0.1 
+0.1 


Total 




+78.4 
-118.4 


+34.2 
-38.6 


+21.4 
-23.2 


+18.3 
-18.1 


+17.9 
-16.6 



Table 5: Central results for ms{Mr) corresponding to the PDG average for \Vus\ extracted 
from the different moments, as well as ranges for the main input parameters and resulting 
uncertainties for rUs- For a detailed description see the discussion in the text. 

The result of our second analysis is displayed in table ^. Like in table ^, again the 
first row contains the individual moment results for central input parameters, then the 
variations of while varying the dominant input parameters are tabulated, and the last 
row gives the combined uncertainty by adding all errors in quadrature. The variation of 
the various input parameters proceeds in complete analogy to the unitarity case discussed 
above. Again performing a weighted average of the individual moment results treating the 
uncertainties like for the unitarity \Vus\, the resulting strange quark mass turns out to be 

ms{Mr) = 107.0 ± 17.9 MeV ^ m,(2GeV) = 103 ±17 MeV. (4.5) 

We observe that for the lower value of \Vus\, we also find to be lowered by roughly one 
sigma, but nevertheless certainly compatible with the result of eq. ( |4.4|) . As table ^ shows, 
here the k dependence of is weaker than for the unitarity case, the reason being that 
now the central lies much closer to the value preferred by the scalar and pseudoscalar 



13 



sum rules. Then the replacement of theoretical by phenomenological spectral functions 
only shifts the resulting by a smaller amount. 



5 A novel route to W, 



us 



Since the quantities 6R'^'' of eq. ( p.lO|) are strongly dependent on \ Vus\, especially for low k, 
now we want to turn the table and investigate the potential of determining IV^^I from the 
SU(3)-breaking r sum rule. For this novel approach, we require a value for the strange mass 
from other sources as an input so that we are in a position to calculate from theory. 
Like in section 3, we have used ms(2 GeV) = 105 ± 20 MeV which is compatible with the 
recent sum rule determinations |^,|y,|21|,^ and also with the lattice results ||41| , ^^ ,^.P] 



Parameter 


Value 


(0,0) 


(1,0) 


(2,0) 


(3,0) 


(4,0) 




central 


0.229 


0.270 


0.320 


0.375 


0.436 


m,(2GeV) 


105 ± 20 MeV 


+0.028 
-0.023 


+0.043 
-0.036 


+0.058 
-0.048 


+0.075 
-0.062 


+0.094 
-0.077 




2x0(0^) 
no 0(03) 


-0.004 
+0.004 


+0.004 
-0.004 


+0.013 
-0.013 


+0.024 
-0.024 


+0.038 
-0.038 




1.5 
0.75 


+0.007 
-0.013 


+0.001 
-0.012 


-0.007 
-0.007 


-0.018 
+0.000 


-0.032 
+0.013 




0.334 ±0.022 


-0.007 
+0.006 


-0.005 
+0.003 


-0.001 
-0.001 


+0.004 
-0.007 


+0.011 
-0.014 


{ss) / (uu) 


0.8 ±0.2 


+0.002 
-0.002 


+0.010 
-0.010 


+0.018 
-0.018 


+0.024 
-0.024 


+0.030 
-0.031 


^scalar 
rus 


see text 


+0.004 
-0.004 


+0.002 
-0.002 


+0.001 
-0.001 


+0.001 
-0.001 


+0.001 
-0.001 


fK 


113 ± 1 MeV 


+0.002 
-0.002 


+0.002 
-0.002 


+0.002 
-0.002 


+0.002 
-0.002 


+0.002 
-0.002 


fx (14.60) 


21.4 ±2.8 MeV 


+0.002 
-0.002 


+0.001 
-0.001 


+0.001 
-0.001 


+0.000 
-0.000 


+0.000 
-0.000 


Total 


error 


+0.030 
-0.028 


+0.044 
-0.040 


+0.062 
-0.054 


+0.083 
-0.074 


+0.107 
-0.098 



Table 6: Central results for extracted from theory, as well as ranges for the main input 
parameters and corresponding resulting uncertainties. 

Our results for the theoretical prediction of SR^'' are presented in table ^. Like our 
results for the strange mass, the first row shows our central values for 6R^'' and below the 
uncertainties resulting from varying the various input parameters within their ranges, now 
including the strange mass, are listed. Again, the last row gives the total uncertainty by 
adding all errors in quadrature. One observes that the errors on 6R'^^ for all moments are 
dominated by the present uncertainty in the strange quark mass, but also the uncertainties 
resulting from higher order QCD corrections play some role. 
^Further references to original works can be found in these reviews. 
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Since experimental results for the individual i?-ratios Rr = 3.642 ± 0.012 and Rr,s = 
0.1625 ± 0.0066 are available, we can employ our result for the (0, 0) moment 



6Rr,th = 0.229 ±0.030, (5.1) 

to calculate the value of | Vus I • Assuming unitarity of the CKM matrix in order to express 
the CKM matrix element \Vud\ in terms of \Vus\, we find 

\Vus\ = 0.2179 ± 0.0044exp ± 0.0009th = 0.2179 ± 0.0045 . (5.2) 

The first given error is the experimental uncertainty due to Rr and Rr.s whereas the second 
error stems from the theoretical value for 6Rr- Our result for \Vus\ is in good agreement 
with the PDG average although at present the uncertainty is roughly twice as large and 
it is only one sigma away from IV^^I as obtained from the CKM unitarity fit. However, 
the current uncertainty is dominated by the experimental uncertainty in Rt,s, and if this 
shrinks by a factor of two, also a corresponding reduction in the uncertainty for \Vus\ would 
be reached. Then the determination of \Vus\ from the r sum rule could reach a precision 
better than the current world average. In addition, if the value of will become better 
known in the future, there is also room for improvement on the theory side. 

A second possibility to calculate \Vus\ from r decays into strange particles with net 
strangeness would be to directly employ the theoretical prediction for Rt.^s/\Vus\- Estimat- 
ing the uncertainties in analogy to table ^, we obtain 

R^^s/\V^^\ = 3.395 ±0.069. (5.3) 

Since for Rr,s^ contrary to 5Rr, also the D = contribution S^^^ is present, in this case 
the uncertainty is dominated by the perturbative QCD error, although also the error 
on the strange quark mass matters somewhat. Together with the experimental value 
Rr^s = 0.1625 ± 0.0066, again we are in a position to extract \Vus\ with the result: 

iKsl = 0.2188 ± 0.0044exp ± 0.0022th = 0.2188 ± 0.0049 . (5.4) 

We observe that our findings of eqs. ( |5.2| ) and ( ^.4| ) are completely compatible. The 
experimental uncertainty in the latter result is identical to the one of ( |5.2| ), whereas the 
theory error, due to the contribution from perturbation theory, is larger. To conclude, 
when improved experimental data will be available, we expect that with our first approach 
a precise determination of \Vus\ from 6Rt- can be achieved. 
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6 Comparison with other analyses 



Our determination of the strange quark mass from hadronic r decays presented above can 
readily be compared to the recent analysis fl^, the most notable difference being our re- 
placement of the badly behaved longitudinal contributions by phenomenological parametri- 
sations. In ref. the final result for the strange mass was rris^MT.) = 120]'^26 assuming the 
value for \Vus\ as obtained from a fit to the unitarity of the CKM matrix. This finding is in 
perfect agreement to our result ( |4.4| ) in this case, although through the usage of the phe- 
nomenological scalar and pseudoscalar spectral functions, a reduction of the uncertainties 
could be achieved in our investigation. 

Other analyses of from hadronic r decays, which have been published in the last 
years, can be found in refs. ||lT],[T^,|l3|,|l4|,[l^,[l6|. A detailed comparison of these analyses 
has recently been worked out by Maltman in ref. |^8|. There it was observed that, once 
all analyses are updated by using common values for R^- s and of the CKM matrix element 
|V^s|, the resulting strange mass values all cluster around ms(2GeV) = 115 MeV for a 
unitarity \ Vus\ and ms{2 GeV) = 102 MeV for the PDG average, being in perfect agreement 
to our findings of eqs. ( [4.4| ) and ( [1.5[ ) respectively. 

Sum rule determinations of the strange quark mass have also been performed on the 



basis of the divergence of the vector or axialvector spectral functions alone ||59|, |60i|21| ,p8 
Being directly proportional to m^, these sum rules are very sensitive to the strange mass 
and the relevant spectral functions are identical to the ones employed to replace the scalar 
and pseudoscalar longitudinal contribution to the r sum rule. The outcomes of these 
analyses were 171^(2 GeV) = 100 ± 12 MeV for the pseudoscalar finite energy sum rule pT| , 
and ms{2 GeV) = 99 ± 16 MeV in the case of the scalar Borel sum rule [^. Again, these 
results are in very good agreement to our value for of eq. ( [4. 51) for the PDG average for 



\Vus\, whereas in the case of the unitarity \ Vus\ the result from the r sum rule is roughly one 
sigma higher. However, at present, due to the large uncertainties no further conclusions 
from this slight discrepancy can be reached. 

The status of the extraction of from the hadronic e~^e~ cross section is less clear. In 
an analysis which took into account flavour-breaking differences, analogous to the r sum 
rule, Narison ^ obtained m,(2GeV) = 144 ±21 MeV. In refs. |3|,|6|] it was then pointec. 
out that large isospin breaking corrections significantly lower the result for the strange mass 
to about ms{2 GeV) = 95 MeV and yield considerably larger uncertainties of the order of 
45 MeV. A recent improved analysis taking into account further SU(3)-breaking differences 
of hadronic correlators now yields a value of ms{2 GeV) = 129 ±24 MeV Nevertheless, 
in our opinion further work in this channel is needed, before a definite conclusion can be 
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Figure 1: Summary of results for ms{2 GeV). 



reached [|63| . 

Recent reviews of determinations of the strange quark mass from lattice QCD have 
been presented in refs. ||4l|, ^ , with the conclusions ms(2GeV) = 108 ± 15 MeV 

and ms(2GeV) = 90 ± 20 MeV in the quenched and unquenched cases respectively. In 
the unquenched case the error in these results is dominated by the uncertainty resulting 
from dynamical fermions, whereas the calculations of in the quenched theory, based for 
example on the kaon mass, are already very precise. However, also in the quenched lattice 
simulations additional uncertainties arise due to the fact that strange masses extracted 



from different quantities display systematic deviations [Q. Nevertheless, the agreement 
between lattice QCD and QCD sum rule determinations of is already rather satisfactory. 

In figure |I], we display a summary of the status of strange mass determinations from 
QCD sum rules and lattice QCD. The two open circles represent the results of this work 
given in eqs. ( [4.4| ) and ( [4.5[ ). The full circles correspond to other sum rule determinations 
discussed in this section, while the full squares are averages for lattice QCD simulations in 
the quenched and unquenched case. Finally, the full diamond corresponds to our average 
value ms{2 GeV) = 105 ± 20 MeV which we have used in sections 3 and 5 as an input. 
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7 Conclusion 



The most severe drawback of determinations of the strange quark mass from the flavour 
SU(3)-breaking difference ( p.lO|) so far has been the bad behaviour of the longitudinal 
contributions ( |3.1| ) within perturbative QCD. This drawback could be circumvented by re- 
placing the badly behaved longitudinal contributions with corresponding phenomenological 
parametrisations of the relevant spectral functions. To gain confidence in this treatment, 
in section 3 we verified that the scalar and pseudoscalar r-like finite energy sum rules are 
satisfied by themselves within the uncertainties, which however are large for the theoretical 
expressions. 

Already in previous analyses it was found that the extracted strange mass depends 
sensitively on the value of the CKM matrix element \Vus\ which is used in eq. ( 2.101) . Thus, 



in our numerical analysis of section 4, we decided to consider two different cases: in the 
first one, like in ref. [0, we employed the value \Vus\ = 0.2225 ± 0.0021 which results from 



a fit to the unitarity of the CKM matrix |^ . At a scale of 2 GeV our resulting strange 
mass (|4.4| ) was found to be 

m,(2GeV) = 117±17MeV. (7.1) 

On the other hand using the PDG central average \Vus\ = 0.2196 ± 0.0026, which is based 
on the analyses |^4|,|5^,^, the strange mass obtained is lowered by roughly one standard 



deviation and takes the value 

m,(2GeV) = 103±17MeV. (7.2) 

One should, however, remark that these results are obtained from an average over the 
different {k, 0) moments which have been analysed and that the individual results for the 
strange mass display a sizeable dependence on the number of the moment. Nevertheless, 
the deviations to our central averages are within one sigma for all individual moments. 

Taking advantage of the strong sensitivity of the flavour-breaking r sum rule on |V^<j|, in 
section 5 we inverted the line of reasoning and determined the CKM matrix element \Vus\ 
from the same sum rule by assuming an average for the strange quark mass as extracted 
from other sources. The result thus obtained is 

\Vus\ = 0.2179 ±0.0045, (7.3) 



where the error is largely dominated by the experimental result Rr^s = 0.1625±0.0066 
A reduction of this uncertainty by a factor of two would therefore result in a corresponding 
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reduction in the error on \Vus\, which would lead to a determination more precise than 
the current PDG average. Furthermore, precise experimental measurements of R^''^ and 
the SU(3)-breaking differences 6R'^^ would open the possibility to determine both and 
\Vus\ simultaneously. This can hopefully be achieved with the BABAR and BELLE r data 
samples in the near future. 

Such an independent access on \Vus\ would be extremely welcome in view of the fact 
that at present the situation about unitarity in the first row of the quark-mixing matrix is 
slightly confusing. Whereas with the PDG average \ Vus\ = 0.2196 ±0.0026, based solely on 
the analysis of Kes decays, unitarity would be violated at the two sigma level, preliminary 
results of the E865 experiment at Brookhaven National Laboratory [0] would correspond 
to \ Vus\ = 0.2278 ± 0.0029 thus being completely compatible with the unitarity of the 
CKM matrix.^] This preliminary value would, however, be about two standard deviations 
away from our result for \Vus\ presented above. In addition, one should keep in mind 
that also shifts in the value for \Vud\ are certainly possible. Anyhow, with the upcoming 
improvements on this subject, we are confident that the question about unitarity in the 
CKM matrix can be resolved in the near future. 
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